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ABSTRACT 

The symmetric peak observed in linear polarization in the core of the solar sodium Di line 
at 5896 A has remained enigmatic since its discovery nearly two decades ago. One reason is 
that the theory of polarized scattering has not been experimentally tested for multi-level atomic 
systems in the relevant parameter domains, although the theory is continually being used for the 
interpretation of astrophysical observations. A laboratory experiment that was set up a decade 
ago to find out whether the Di enigma is a problem of solar physics or quantum physics revealed 
that the Di system has a rich polarization structure in situations where standard scattering 
theory predicts zero polarization, even when optical pumping of the m state populations of the 
hyperfine-split ground state is accounted for. Here we show that the laboratory results can be 
modeled in great quantitative detail if the theory is extended to include the coherences in both 
the initial and final states of the scattering process. Radiative couplings between the allowed 
dipole transitions generate coherences in the initial state. Corresponding coherences in the final 
state are then demanded by a phase closure selection rule. The experimental results for the well 
understood D 2 line are used to constrain the two free parameters of the experiment, collision rate 
and optical depth, to suppress the need for free parameters when fitting the Di results. 

Subject headings: line: profiles - methods: laboratory: atomic - polarization - radiation mechanisms: 
general - scattering - Sun: atmosphere 


1. Introduction 


The development of new theories will often go 
astray unless guided and tested by observations 
or experiments. In the early stages of quan¬ 
tum physics experiments with polarized scatter¬ 
ing helped guide the development of the theory. 
The experiments by Wilhelm Hanle in Gottingen 
( Hanlel 1924) led not only to the discovery of 
the well used Hanle effect but at the same time 
demonstrated the concept of coherent superposi¬ 
tion of quantum states and the partial decoher¬ 
ence caused by external magnetic fields. The new 
quantum theory was built on these concepts. 


The literature on experiments with the scat¬ 
tering of polarized light came to an abrupt end 
around 1935, apparently because the topic was 


considered exhausted and less rewarding than 
other areas, and because the observed phenomena 
were believed to be sufficiently understood. The 
theoretical understanding covered the parameter 
domain accessible at that time, which in terms of 
polarimetric sensitivity was unbelievably crude in 
comparison with what is now possible with mod¬ 
ern technology. Therefore the quantum scattering 
theory never got experimentally tested over a wide 
range of weak polarization phenomena that are ac¬ 
cessible with current instrumentation. The theory 
is nevertheless being systematically applied for the 
interpretation of astrophysical observations in this 
untested parameter domain. 


With the implementation in the 1990s of 
a new techn ology (ZIMPOL) in sol ar spectro- 
polarimetry ( Povell 19951 2001 : Gandorfer et al 
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2004T) . which allowed the recording of the Sun’s 
spectrum with a polarimetric precision of 10 5 , 
the linearly polarized spectrum (Stokes Q/I) ob¬ 
served near the solar limb was found to be as 
spectrally structured as the intensity spectrum, 
but the spectral structures were totally different, 
it wa s like uncovering a new spe ctral face of the 
Sun ( Stenflo fc Keller 19961 1 1997 ). It is therefore 
referred to with the term ‘Second Solar Spectrum’, 
or with the acronym SS2 (implying that the or¬ 
dinary intensity spectrum is SSI). The origin of 
this polarized spectrum is coherent scattering pro¬ 
cesses in the Sun’s atmosphere. Magnetic fields 
do not produce this spectrum, but they modify it 

SS2 is a playground for the Hanle effect. 

Through intense theoretical efforts most of 
these unfamiliar structures could be identified 
and understood within the framework of standard 
quantum scattering theory, revealing remarkable 
polarization signatures of quantum interference 
between states of different total angula momenta 
(J quantum numbers), lryperfine structure and 
isotope effects, optical pumping, molecular scat¬ 
tering, and partial frequency redistr ibution effects 
in polarized radiative transfer (cf. IStenfio 2004 . 
for a review). 

A culmination of the theoretical work to build 
a comprehensive quantum-mechanical foundation 
for the radiation-matter interactions occurred 
with the publi cati on of the monu menta l mono¬ 
graph by Landi Degl’Innocenti fc Landolfil ( 2004 ). 
A regular series of International Workshops has 
been devoted to the subject of ‘Solar Polariza¬ 
tion’ since 1995, the late st, no. 7, in Kunminj 


China, September 2013 (INasrendra et al.1 1201 


The physics of polarized scattering has been a 
central topic of these Workshops. 

In spite of these various efforts, one kind of ob¬ 
served feature in the Second Solar Spectrum has 
remained enigmatic and stubbornly resisted all 
attempts to explain it, namely the weak linear po¬ 
larization peak observed at the center of the Nal 
Di line at 5896A an d the Ban Di line at 4934A 
( Stenflo et al.ll2000al ibl) . although attention to the 
enigmatic sodium Di feat ure was drawn already i n 
the 1996 paper in Nature ( Stenflo fc Keller 1996h . 
The problem is that Di type lines, which have 
angular momentum quantum numbers J = 1/2 
for both the lower and upper atomic level is ex¬ 
pected to be intrinsically unpolarizable, represent¬ 


ing polarization ‘null’ lines. However, since both 
sodium and barium have nuclear spin 3/2, both 
the lower and upper levels are split into hyperfine 
structure states with total angular momentum 
quantum numbers F = 1 and 2. A number of 
attempts have been made to explain the observed 
Di polarization in terms of optical pumping of 
the hyperfine structur e level s of the ground state 


j LandCDegTInnocenti 1998|; Casini et al. 2002; 


Casini fc Manso Sainzj 20051) . but all of them have 
failed because the predicted Di polarization is too 
small by about two orders of magnitude and has 
the wrong symmetry (anti-symmetric polarization 
profile i nstead of the obser ved symmetric one) 
( Trujillo Bueno et al. 20021: Kerkeni & Bommicr 


20021: iKlement fc Stenfloll2003i)~ 

To find out if this enigma is a problem of solar 
physics or of quantum physics, a laboratory exper¬ 
iment was set up to explore the polarization struc¬ 
ture of Di scattering under controlled conditions. 
The experiment unambiguously showed that there 
is indeed a problem with the quantum theory of 
polarized scattering, since the theory predicts null 
results when the e xperiment revealed a r i ch po lar- 
ization structure ( Tlialmann et al.l 20061 20091) . 

In search for ways in which the theory of quan¬ 
tum scattering would need to be extended to ex¬ 
plain what was seen in the laboratory experiment, 
it was realized that if coherences in both the ini¬ 
tial and final ground states could be included in 
the theory for the scattering process, then the 
atomic system would possess a much richer reso¬ 
nant structure that would be able to gene rate po¬ 
lariza tion signatures of the observed kind ( Stenflo 
20091) . 

In the present paper we develop these ideas into 
a phenomenological extension of current scatter¬ 
ing theory and use it to successfully model the 
observed Di polarization in great quantitative de¬ 
tail. Through use of experimental data in the well- 
understood D 2 line to constrain the values of the 
collisional and optical depth parameters y c and r 
it becomes possible to avoid the use of adjustable 
parameters when fitting the Di data. 

2. Laboratory experiment 

After initial attempts failed to do the labora¬ 
tory exploration of Di scattering physics for the 
sodium Di line with the use of a broad-band light 
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source, we realized that a tunable laser was re¬ 
quired to achieve the needed S/N ratio with clean 
separation between the D 2 and Di transitions and 
insignificant stray light. In addition the tuning 
gives us the precise shapes of the polarization pro¬ 
files. 


We then chose to do the experiment for potas¬ 
sium instead, since solid-state tunable lasers are 
not available for the sodium wavelength range. 
The Kl D 2 7665 A and the Kl Di 7699 A lines 
have the same quantum number structure as the 
Nai D 2 5890 A and the Nai Di 5896 A lines, in¬ 
cluding nuclear spin and hyperfine structure, so 
it is the same physics that is being tested with 
potassium as for sodium. 


The set-up has been de scribed in detail 
( 2006 12009.; together with 


Thalmann et al. 


by 

exam¬ 
ples of the experimental results. Here we provide a 
summary of what is relevant for the understanding 
of the theoretical modeling of the data. The heart 
of the system is a cross-shaped glass cell that con¬ 
tains the potassium gas. The cell was tailor-made 
for our experiment by the late Alessandro Cac- 
ciani in Rome, who had perfected the art of build¬ 
ing such cells for use in magneto-optical narrow- 
band filters used in helioseismology and solar mag- 
netometry ( Cacciani fe Fofi 1978t Cacciani et al 


1997]). Metallic potassium in the stem of the cell 


is heated to give potassium gas at a temperature 
of 100° C at the cell center, where the scatter¬ 
ing takes place. An argon buffer gas is used to 
suppress diffusion of the potassium vapor towards 
the cooler entrance and exit windows, to avoid de¬ 
posits that could make the windows opaque. 


The spectral band width of the laser beam is 
not exactly known, but it is at least less than 
2 mA. This is much smaller than the separation 
between the F = 1 and F = 2 hyperfine levels of 
the ground state, which is 9.1mA, but it is com¬ 
parable to the corresponding hyperfine splitting of 
the upper state of the K Di line, which is 1.14 mA. 


The expanded laser beam is passed through po¬ 
larizers that allow us to select one of six states of 
100% fractional polarizations in Stokes Q, U, or 
V before the beam enters the cell: 100% ±Q/I , 
±U/I, or ±V/I. Positive Stokes Q is defined to be 
linear polarization oriented perpendicular to the 
scattering plane. The scanning can be done with 
either of two separate laser heads, one for tuning 
around the K D 2 resonance, the other for tuning 


around the K D[ resonance. 

In the output arm that receives light scattered 
by the cell at 90°, a piezoelastic modulator fol¬ 
lowed by a linear polarizer converts the polariza¬ 
tion information (Stokes Q, U , or V) into intensity 
modulation, at 84 kHz for Q and U , at 42 kHz for 
V, while Stokes I is represented by the unmodu¬ 
lated (DC) signal. A photomultiplier feeds the sig¬ 
nal to a lock-in amplifier, which demodulates the 
AC component. The combination of 6 alternative 
input states with three alternative output states 
for the polarization gives us 18 possible combina¬ 
tions. For each of them the laser tuning gives us 
the precise profile shape with mA resolution. 

Furthermore a set of Helmholtz coils mounted 
on the cell allows us to impose on the cell cen¬ 
ter an external magnetic field with any strength 
in the range between ±30 G, with an orientation 
that can be chosen to be in either of the three 
spatial directions: perpendicular to the scattering 
plane, parallel to the input beam, or parallel to 
the scattered beam. 

Because the available parameter space of the 
experiment is so large, time and manpower only 
allowed the experiment to be carried out for some 
strategically selected combinations among the var¬ 
ious possibilities. In the present paper we focus on 
the interpretation of the transverse field case (field 
perpendicular to the scattering plane) for the two 
combinations (1) Input 100% Q/I polarization, 
output detection of Stokes Q , and (2) input 100 % 
V/I polarization, output detection of Stokes V. 
This selection was done not only because these 
cases were most thoroughly measured, but because 
conventional scattering theory predicts null results 
for case (1) and nearly null for case (2). This leads 
us to identify the observed polarization effects for 
these cases as due to previously overlooked, ne¬ 
glected physics. 

3. Physics of polarized scattering 
3.1. Theoretical overview 

Let us consider a scattering transition between 
the magnetic substates labeled a (initial state), 
b (intermediate), and / (final) and introduce the 
area-normalized complex profile function for ra- 
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diative absorption: 


= 


2/i 

Uba - u ~ *7/2 ’ 


(1) 


where 

Uba = (E b - E a )/h (2) 

represents the transition frequency between the 
energy levels of the upper and lower magnetic sub¬ 
states, while ui is the frequency of the incident ra¬ 
diation. Then the Kramers-Heisenberg formula for 
the probability amplitude for the scattering pro¬ 
cess a —>• b —> f may be written as 


w a p ~(/|r-e a |fe)(&|r-e /3 |a) <F m/ia (3) 


(cf.|St en flol ll994L IToOSh . 


The matrix transition elements contain the 
scalar products between the dipole moment op¬ 
erator (proportional to the position vector r) and 
the unit vectors ep and e a of the incident and 
scattered radiation field, respectively. Vector r 
is decomposed in the three spherical vector com¬ 
ponents labeled by q = 0, ±1, which equal the 
difference Am between the magnetic quantum 
numbers of the lower and upper states (with the 
quantization axis defined to be along the magnetic 
field vector), while the polarization components of 
the radiation field are real linear unit vectors that 
reside in the plane that is perpendicular to the 
propagation direction. 


For classical scattering the probability ampli¬ 
tude is represented by the Jones matrix 


w, 


o r^j (I) 

yp ^ — q c q : 


(4) 


where the e factors embody th e scal ar products 
between r and e a> p (cf. ISten flol 11994 . With our 
definition of q , we have <f>_ g = & ba for the case 
of a J = 0 —y 1 absorption transition. In the 
general quantum-mechanical case we can write 
the Kramers-Heisenberg amplitude in the follow¬ 
ing analogous form, to bring out the similarlities 
and differences with respect to the classical case 
as transparently as possible: 


'Waft ^ *&ba &q> £q tfb tba ■ (b) 

Here t ba and tf b represent the transition ampli¬ 
tudes for absorption and emission between the re¬ 
spective magnetic substates. They are all real 


quantities, but their relative signs are of impor¬ 
tance. All the other factors in Eq. d5j) are in gen¬ 
eral complex-valued quantities. 

Since the magnetic quantum numbers for the a 
and / states may differ, the corresponding q in¬ 
dices may also differ, which is the reason for using 
a prime for q in the case of the emission transition. 
In the classical case q' = q, and the three reso¬ 
nances (corresponding to the three q) have equal 
probability amplitude and are driven in phase by 
the incident radiation field. In the quantum case 
the scattering amplitudes get weighted by the rel¬ 
ative complex amplitudes c 0 of the initial state a. 

To go from the scattering probability ampli¬ 
tudes w a fi , which represent the elements of a com¬ 
plex 2x2 matrix w, to obtain the scattering prob¬ 
abilities for the different polarization states, one 
forms the bilinear products w a p w* a ,p ,, which are 
elements of a tensor product of type w (g> w* with 
appropriate summations over the various atomic 
states. 


The standard approach in quantum scattering 
theory has been the ‘sum over histories’: To ob¬ 
tain the total probability amplitude for scattering 
from state a to state /, one sums over the proba¬ 
bility amplitudes for all the possible intermediate 
states, while keeping the initial and final states a 
and / fixed. One then forms the bilinear product 
to obtain the probabilities, and afterwards sums 
over all the possible initial and final states a and 
/. This results in a complex 4x4 coherency matrix 
W of the form 


W = 5Il c “| 2 51 (55 w ) 0 (55 w *) 


(IStenflolll99lll998h . 


( 6 ) 


There is a very fundamental distinction be¬ 
tween the various summations here: While the 
sum over the intermediate states is coherent (done 
as a linear superposition of the probability ampli¬ 
tudes), the sums over the initial and final states 
are incoherent (done over the bilinear products 
that represent probabilities), weighted by the rel¬ 
ative populations |c 0 | 2 of the initial states. The 
quantum interference phenomena are then exclu¬ 
sively due to the cross products that arise when 
the sums over b and b r are multiplied with each 
other. 


In the present paper we will present experi¬ 
mental and theoretical evidence that demonstrates 
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that Eq. f6]) with its separation between coherent 
and incoherent summations is incorrect, and that 
instead all summations shoul d be done coherently 
(as previously suggested in lSten fiol l2009ll . The cor¬ 
rect expression should therefore be 

w = (5Z c “ w ) ® c “' w *) • 

abf a'b' f 

Note that it is absolutely essential to make the 
summation coherent not only over the initial states 
(accounting for the non-diagonal density matrix 
elements c a c*,), but also use coherent summation 
over the final states. This rule is demanded both 
by a theoretical selection rule (cf. Sect. 13.41) and 
by the experimental data. If it is obeyed, then we 
can reproduce in great detail the results of the lab¬ 
oratory experiment, while if it is disobeyed we get 
zero polarization in the corresponding test cases. 
The selection rule in Sect. [33] implies that there is 
a strict coupling between the initial and final state 
coherences. 

Since in this case each sum contains a much 
larger number of terms, the physical system will 
contain a much richer set of coherences and quan¬ 
tum interference phenomena, which are prevented 
in the earlier approach to quantum scattering as 
represented by Eq. ([6]). The new terms allow us 
to theoretically reproduce the results of our lab¬ 
oratory experiment in great quantitative detail in 
situations when the old formalism predicts zero 
polarization. 

The practical implementation of the formally 
simple and symmetric Eq. 0 is however difficult, 
both conceptually and technically, since we need a 
way to determine the phase relations between the 
different ground level amplitudes c a . In the ab¬ 
sence of a radiation field the relative phases are 
random. In the presence of an incident radia¬ 
tion field the transverse oscillations of the electric 
vector drive the dipole resonators of the atomic 
system. It is this synchronized driving that pro¬ 
duces phase locking between combinations of c a 
level amplitudes. The heart of the problem is 
how to deal with this phase locking, which must 
be distinguished from conventional optical pump¬ 
ing, although it may superficially resemble it. To 
conceptually understand the fundamental distinc¬ 
tion between phase locking and optical pumping 
we need a non-traditional way to look at the na¬ 
ture of quantum transitions, an issue that will be 


addressed in the next subsection. 

Before proceeding with the theoretical issues, 
let us already at this stage present experimental 
evidence that backs up our previous statements. 
In Fig. |T]we show the case of 90° scattering of K 
Di radiation when the incident light is 100 % lin¬ 
early polarized perpendicular to the plane of scat¬ 
tering (incident Stokes Q/I = 1), and we measure 
Stokes Q of the scattered radiation. The laser has 
been tuned across the Di resonance, which gives 
us the profile shape. Standard scattering theory 
(Eq. J6])) predicts exactly zero Q polarization for 
all wavelengths, while with Eq. m we are able to 
reproduce the observations in great quantitative 
detail. Note that the observed Q is negative, i.e., 
the scattered radiation is weakliy linearly polar¬ 
ized parallel to the scattering plane. 

In the following sections we will explain what 
explicitly goes into these calculations, how the the¬ 
ory is extended to include the new physical effects, 
and how observations of the well-behaved D 2 sys¬ 
tem are used to constrain the theory and enforce 
consistency. An objective is to provide convincing 
evidence that quantum scattering has a far richer 
coherency structure than allowed by the old for¬ 
malism of Eq. ©, and to indicate the direction 
in which the theory should be extended to cover 
these phenomena. 

To conclude this subsection we give the well- 
known expression for how to obtain the Mueller 
matrix M that describes scattering of the Stokes 
vector from the complex coherency matrix W: 

M = T WT -1 . (8) 


The complex 4x4 matrices T, T -1 in Eq.© 
are purely mathematical tran sformat ion matrices 
without physical contents (cf. jstenflo 1994 ). 


3.2. Wave packet interaction: Quantum 
blip or radiation bath ? 

An obstacle to a more complete understanding 
of the interaction of the atomic system with the ra¬ 
diation field has been the often prevailing miscon¬ 
ception that atomic transitions are instantaneous 
events, quantum ‘blips’. To make our conceptual 
discussion of the interaction process as transparent 
as possible we will refrain from applying the sec¬ 
ond quantization formalism with the quantization 
of the electromagnetic field and instead use a semi- 
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Fig. 1.— Polarized Di scattering by potassium 
when the incident radiation is 100 % linearly polar¬ 
ized perpendicular to the scattering plane (Q/I = 
1). The radiation scattered at 90° is then partially 
polarized parallel to the scattering plane (output 
Q negative). The observations (dashed curves) 
can be reproduced by the theory (solid curves) 
when the coherences in both the initial and fi¬ 
nal states are accounted for as in Eq. (0. If the 
coherences in the final state would be neglected, 
we would get zero polarization for all wavelengths, 
regardless of how unbalanced the m state popula¬ 
tions are. 


classical approach in which photons are treated as 
exponentially damped wave packets. 


The misconception of instantaneous quantum 
jumps has been amplified by the circumstance 
that a theory for polarized optical pumping is 
only available under the assumption that the 
incident radiation is broadband in frequency 
(the fl at-spectrum app roximation, cf. p. 257 in 
Lancli Degl’Innocenti & Landolfi 20041) . Although 


the flat-spectrum assumption only requires the 
radiation to be spectrally flat over ranges compa¬ 
rable to the Doppler width of spectral lines, let us 
for conceptual clarity first consider the extreme 
case of complete flatness. With Heisenberg’s un¬ 
certainty principle AEAt ss h such a case implies 
that At = 0, i.e., the impinging wave packets 
have zero coherence length. Such wave packets do 
not exist in nature. Furthermore, this limitation 
contradicts the principle that each wave packet 
can be represented as a linear superposition of its 
monochromatic Fourier components. 


While the flat-spectrum approximation does 
not imply this extreme case of zero coherence 
length, it does imply wave packets with a coher¬ 
ence length that is about two orders of magni¬ 
tude smaller than actual wave packets, since the 
Doppler width is typically two orders of magni¬ 
tude larger than the radiative damping width. For 
this reason the implications of this approxima¬ 
tion are qualitatively of the kind indicated by our 
discussion of completely flat radiation. Thus, in 
terms of hypothetical (non-existing) broad-band 
wave packets the pumping scenario is a sequence 
of almost discrete pumping events, in which each 
quantum jump due to radiative excitation is nearly 
instantaneous. During the relatively long waiting 
times between each discrete event the atomic sys¬ 
tem is subject to the depolarization effects of col¬ 
lisions and magnetic field s (cf . the more extended 
discussion of this topic in Stenflol 2011 ). 


It is only meaningful to talk about quantum 
‘blips’ when a wave packet is recorded in a detector 
(collapse of the wave function). Before this hap¬ 
pens the wave packets are not small entities but 
huge, with coherence lengths of several meters (in 
the absence of collisional effects). The radiative 
damping constant of potassium Di and D 2 is ap¬ 
proximately 7 = 3.8 x 10 7 s -1 (determined by the 
inverse life time of the excited state in the colli¬ 
sionless case), implying that the damped emission 
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process takes place over a time of order I/ 7 , which 
is 10 7 times longer than the time scale for the oscil¬ 
lation of the electric vector. During this long time 
light travels approximately 8 m, corresponding to 
10 7 wavelengths. 


With this enormous separation of time scales 
wave packets represent (from the point of view 
of the atomic system) a radiation bath. During 
the interaction with the atomic system the various 
dipole resonances in the atom are driven 10 ' times 
by the electric field of the single wave packet. This 
gives plenty of time to establish synchronization 
between the phases of the various dipole oscillators 
represented by the potential radiative transitions 
between lower and upper magnetic substates. In 
the case of the Di atomic system there are 36 dif¬ 
ferent allowed combinations of dipole transitions 
between lower and upper m states. Since all of 
them are subject to the driving force from the 
same external oscillating electric field (provided 
by a single incident wave packet), phase synchro¬ 
nizations will be established, depending on the po¬ 
larization of the incident wave packet. 


Both the lower and upper J = 1/2 states of 
the Di atomic system is split due to the nuclear 
spin into total angular momentum states F = 1 
(with three magnetic substates) and F = 2 (with 
five magnetic substates). There are thus 8 mag¬ 
netic substates in each of the lower and upper 
state. Various combinations of the 36 allowed 
dipole transitions have common m states (either 
lower or upper ones). For instance, a Am = 0 
transition from upper state F = 1, m = 0 couples 
to both F = 1, m = 0 and F = 2, m = 0 in the 
lower state. Since both the upwards and down¬ 
wards transitions are jointly driven by the same 
oscillating electric field, phase synchronization be¬ 
tween the two F = 1, ?n = 0 and F = 2, m = 0 
lower states will be produced, which then leads to 
new polarization effects through Eq. ©■ 


Let us for completeness mention here that the 
fully quantized formu lation of polarized scattering 
that was pioneered bv lBommler ( 1997al lbli does not 
make use of the flat-spectrum approximation. Al¬ 
though the method is general with the potential 
of being extended to the multi-level case, it has so 
far only been implemented for the two-level case. 


3.3. Coherences of the atomic system 

We have in the previous subsection seen that 
the wave packet interaction is not associated 
with instantaneous quantum jumps, but that the 
atomic system with its various dipole resonators 
(36 in the Di case) is immersed in a radiation 
bath (provided by each single wave packet), which 
drives the atomic oscillations and couples combi¬ 
nations of oscillators. 

This viewpoint is largely prohibited when the 
flat-spectrum approximation is used, as indicated 
in the previous section. There are nevertheless 
good reasons why this approximation has repre¬ 
sented the established approach, because by re¬ 
moving it we open a ‘Pandora’s box’ of coherency 
phenomena that we do not know how to properly 
deal with. However, by keeping the lid shut on 
this Pandora’s box we prevent access to the world 
of coherency phenomena in which a solution of the 
Di enigma is to be found. 

In the present paper we open the lid to examine 
the phenomena in this new world and try to show 
how they are responsible for the observed polar¬ 
ization features in the laboratory experiment. In 
view of the complexity of the interconnected co¬ 
herences our initial approach has to be heuristic 
and phenomenological in nature rather than rep¬ 
resenting an attempt at a formulation of a new 
general theory. Still we go to great lengths to 
make the phenomenological extension of standard 
scattering theory well defined, internally consis¬ 
tent, based on the soundest possible arguments, 
and avoiding free parameters as far as possible. 
This effort is only possible with the guidance of 
experimental data. The success of our extension 
is that it is capable of modeling the key experi¬ 
mental data in great quantitative detail. 

In the present subsection we will examine the 
nature of the various new coherences. For a con¬ 
ceptual understanding we begin by examining how 
phase synchronization or phase locking plays out 
in classical physics. 

3.3.1. Coherences in classical scattering 

The equation for a classical oscillator in an ex¬ 
ternal magnetic field can be written as 

- ( 2 <Aw l - + Wq r q = ~—E q (9) 
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in component form, when instead of the three 
Cartesian spatial coordinates the three complex 
spherical vector components, labeled with q = 
0, ±1, are used. The q = 0 component is aligned 
along the magnetic field. r q is the component 
of the position vector of the oscillator, ojl = 
eB/(2m) is the Larmor frequency. 7 represents 
the damping due to the radiative energy loss of 
the oscillating system, an d u >p is t he resonance fre¬ 
quency of the system (cf. IStenflolll994l) . 


The advantage of using the complex spheri¬ 
cal vectors rather than real linear vectors is that 
the component equations then decouple and repre¬ 
sent three independent oscillators. Still, although 
the three equations are mathematically uncoupled 
(the equation system is diagonalized), the three 
oscillators are phase locked because the three com¬ 
ponents E q of the driving oscillating electric field 
of the external radiation bath are mutually phase 
locked (since they are components of the same os¬ 
cillating vector). 

The solution of the stationary equation, when 
the oscillator is immersed in an eternal, monochro¬ 
matic radiation bath, leads to the profile function 
which is proportional to the complex refrac¬ 
tive index of the medium. As the dipole radiation 
emitted by these oscillators represents the scat¬ 
tered light, we readily obtain the Jones scattering 
matrix in the form of Eq. ©• 


The frequency redistribution, i.e., the relation 
between the incident and scattered frequencies, 
can be obtained by solving Eq . (|9l) as a time- 
dependent problem ( Bommier fc Stenflol 1999). 
The full solution of this problem is the sum of the 
stationary solution and a transitory solution. The 
collisional effects are generally described in two 
ways: As phenomenologically introduced branch¬ 
ing ratios between the stationary and transitory 
solutions, and as phase truncations of the oscilla¬ 
tors, which lead to depolarization. Here we will 
not deal with the details of such redistribution ef¬ 
fects, since they are not crucial for the key issues 
in the extension of scattering theory. 


However, one aspect that is fundamental to this 
extension and which can be correctly described by 
the classical theory is the partial decoherence of 
jointly driven oscillators that have different reso¬ 
nant frequencies. In the absence of magnetic fields, 
the three classical oscillator components have the 
same resonant frequencies and therefore remain in 


phase for the duration of the radiation bath, there 
is no decoherence. When however there is a fre¬ 
quency split induced by the Larmor frequency, the 
oscillators get progressively out of phase, depend¬ 
ing on the ratio between the time scales of Lar¬ 
mor precession and radiative decay. This partial 
decoherence leads to the polarization effects in the 
scattered radiation that are covered by the term 
‘Hanle effect’. 


The partial decoherence between jointly driven 
resonators, i.e., between jointly driven atomic 
dipole transitions, is by no means limited to deco¬ 
herence between the magnetically split transitions 
(Hanle effect), but applies to the relative phase 
relations and decoherence for any split but jointly 
driven transitions. An example is the dramatic 
polarization effects seen on the Sun due to scat¬ 
teri ng at the Can H and K line system ( Stenflol 
198f)h . two lines of the Di - D 2 type (although 


without hyperfine structure). In this particular 
case the line separation is due to the fine struc¬ 
ture splitting (different J quantum numbers). In 
the potassium case that we are dealing with the 
splitting is due to hyperfine structure (different 
F quantum numbers), and for non-zero magnetic 
fields due to Zeeman splitting as well. 


When forming the bilinear products w a p w* a ,o, 
between the elements of the Jones matrix for the 
scattering probability amplitudes we get interfer¬ 
ence terms, in which the partial decoherence is 
contained in the products ‘hba^'a' between the 
profile factors. An elegant mathematical trans¬ 
formation is to convert such profile products into 
sums: 


QbaH'a' ~ COSC* e la i ($ fea + $£„,) (10) 

( Stenflol 1994. 1998), where the angle a is given by 


tana = 


bJba Wb'a' 
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(ii) 


In the case of magnetically induced splitting a is 
referred to as the ‘Hanle angle’. The expression 
is however quite general and is valid regardless of 
the physical origin of the splitting. 

With the conversion of Eq. m the interfer¬ 
ence effects due to the profile product gets fac¬ 
torized into a frequency-independent Hanle part 
cos a e la and a frequency-dependent part, which 
becomes unity if we integrate over all frequencies 























(since each profile function is area normalized). 
The Hanle part consists of an amplitude factor, 
which in the scattering polarization corresponds 
to depolarization, and a phase factor that corre¬ 
sponds to rotation of the plane of polarization. 

The classical oscillator case corresponds to the 
quantum-mechanical case for a J = 0 -f 1 -> 0 
scattering transition. It also consists of three res¬ 
onators, which link the three m state of the upper 
level to the common, single lower state. In such a 
situation a' = a, so there can be no lower-level co¬ 
herences. In the general case when there are more 
than one lower state, ground-state coherences are 
possible, but since the coherency terms not only 
contain the profile products of Eq. (usd but also 
the level amplitude products c a c*,, these coher¬ 
ences will only play a role if there are correlations 
between the phases of certain combinations of level 
amplitudes. In the absence of an external radi¬ 
ation field the phases are uncorrelated, because 
there is no coupling between the levels. Remain¬ 
ing lower-level coherences will be erased if the time 
between collisions is much shorter than the ‘wait¬ 
ing time’ between succssive radiative excitations. 
The needed synchronization between the phases 
can however be provided by the oscillating elec¬ 
tric field of the wave packet radiation bath, which 
is a source of coupling or phase locking between 
combinations of ground states. 

3.3.2. Coherences in the initial and final states 

The coherences in the initial state are repre¬ 
sented by the off-diagonal elements c a c*, of the 
density matrix. In the absence of outside distur¬ 
bances (external radiation field or collisions), each 
level amplitude oscillates like exp {—iE a t/h). If 
the energy levels E a and E a i are not identical, 
the oscillations will rapidly get out of phase. Let 
us now assume that the phases get synchronized 
by the driving oscillating electric field of an inci¬ 
dent wave packet, but that the initial phase syn¬ 
chronization decays with the damping rate 7 be¬ 
cause the driving electromagnetic wave packet is 
exponentially damped according to this rate. In 
the time domain the density matrix element is 
then exp(— t — r yt). The time derivative of 
this expression must for stationarity be balanced 
by the radiative absorption and emission rates, 
which represent the interaction between the driv¬ 
ing oscillating electric field and the atomic system. 


The balance equation then gives us the expression 
1 /(/ Woo' + 7 ) for the corresponding density matrix 
element, which may more elegantly be written as 


1 

1 + i tan /3 


cos/3 e 1/3 


( 12 ) 


(cf. Eq. (8.89) in Stenflol 1994 ). The angle /3 is 
given by 


tan/3 = 


(13) 


Comparison with Eqs. m and m shows that 
this initial-state coherence is of the same general 
form as ’Hanle type’ coherences between resonant 
transitions. 


If the oscillating electric field of the incident 
radiation has a polarization that leads to radia¬ 
tive coupling between the a and a' states (via ap¬ 
propriate excited states), then we expect the con¬ 
stant of proportionality in Eq. m to be close to 
I c 0 1 | Co'| (the product between the square roots of 
the sublevel populations), since when the phases 
of the level amplitudes are synchronized in the ra¬ 
diation bath, the only source of decoherence is the 
difference in oscillation frequency as expressed by 
the Hanle-type factor in Eq. m- In the absence 
of any radiative coupling the proportionality con¬ 
stant would be zero, with the consequence that 
the respective initial-state coherence does not con¬ 
tribute to the scattering process. In the next sub¬ 
section (Sect. 13.3.31) we will discuss the nature of 
this radiative coupling in some more detail using 
our laboratory experiment for potassium Di - D 2 
as an example. 

As the energy separation between the a and a' 
levels increases, the initial-state coherence terms 
get suppressed towards zero by a ‘depolarization 
factor’ cos 2 /3 represented by the real part of the 
Eq. (fl2l) expression. Because of the symmetry 
between the radiative excitation and deexcitation 
process, we expect there to be a similar suppres¬ 
sion factor governing the final-state coherences. 
There is however a difference between the two legs 
of the a —> b —> f scattering process: the a — > b leg 
is associated with the complex $f, a profile function 
of Eq. (JTJ while the b f leg is not. The imagi¬ 
nary part of Eq. m with its associated minus sign 
is related to the role of the complex profile func¬ 
tion in the oscillation between the a and b states, 
which establishes the radiative coupling. No such 
complex-valued term is involved in the b —> f leg 
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(since the transition amplitudes tfb are real num¬ 
bers). Therefore it is most reasonable to expect 
that the final-state coherences should be governed 
by the symmetric version of Eq. m, namely the 
sum of this expression with its complex conjugate 
(normalized to unity for zero splitting). This nor¬ 
malized sum equals the cos 2 /3 suppression factor, 
with the difference that now /3 is not given by 
Eq. ([bill , but by the same equation if we replace 
aw by uiff>. 

The only valid benchmark that we can use to 
test the validity of these rather heuristic argu¬ 
ments is experimental data. Fortunately the scat¬ 
tering polarization turns out to be extremely sen¬ 
sitive to the choices we make for the ground-state 
coherences. If we omit them, we get zero polariza¬ 
tion, regardless of the relative initial m state pop¬ 
ulations. If we use Eq. G2D but set the final-state 
coherences to be unity (no suppression factor), we 
get much too large polarization effects. If we keep 
an imaginary part of the final-state coherences like 
for the initial state, then we also cannot achieve 
a fit with the data. The same happens if we omit 
the imaginary term in Eq. m or change its sign. 
Also we must only use the correct combinations 
of radiative couplings (cf. Sect. I3.3.U1) . otherwise 
we get results that have no resemblance with the 
observations. Only with the previously described 
choices of expressions for the coherence terms we 
succeed in obtaining a fit at all, and on top of that, 
the fit agrees with the observational data in great 
quantitative detail in a way that is consistent be¬ 
tween Di and D 2 and between the Stokes /, Q, 
and V parameters, including their profile shapes. 

Let us note here that for a given initial-state 
(a, a') combination only certain (/, /') combina¬ 
tions for the final-state coherences are allowed. 
The nature of these selection rules will be ad¬ 
dressed in Sect. 13.41 Without them we find no 
agreement with the observational data. 

3.3.3. Radiative coupling in the D\ - D 2 experi¬ 
ment 

Our analysis here of the laboratory scattering 
experiment will focus on the cases when ( 1 ) the in¬ 
cident radiation is 100 % linearly polarized perpen¬ 
dicular to the scattering plane, while the Stokes I 
and Q parameters of the 90° scattered radiation 
are recorded, and ( 2 ) when the incident radiation 
is 100% circularly polarized, while the Stokes I 


and V parameters of the scattered radiation are 
recorded. The imposed external magnetic field is 
chosen to be vertical (perpendicular to the scat¬ 
tering plane). This direction represents the quan¬ 
tization direction in our theoretical discussion. 

Let us first turn our attention to case (1). Then 
among all the 36 resonators in the Di system 
only those with m a — mb = q = 0 get driven by 
the oscillations of the electric field of the incident 
wave packet, with the consequence that there can 
only be radiative coupling between initial a and a' 
states for which the m quantum number is identi¬ 
cal ( m a i = m a ). The only possibility for two such 
states to represent separate states is to belong to 
different F states. 

Let us as a concrete example see how the two 
states (F = 1, m = —1) and (F = 2, m = —1) are 
radiatively coupled. Each of them couple by ra¬ 
diative excitation to each of the two corresponding 
states of the upper level. Similarly each of the two 
upper states are coupled to the corresponding two 
lower states by the emission process. The oscil¬ 
lating electric field drives the upwards and down¬ 
wards oscillating transitions with equal transition 
strength. 

All the four m = — 1 states (two lower and two 
upper) are therefore mutually phase locked as a 
consequence of the common driving external elec¬ 
tric field. For the same reason that we have phase 
synchronization between the three classical oscil¬ 
lators (which leads to the Hanle factor of Eq. (fTOl) 
for the upper states of a J = 0 —> 1 transition), 
the oscillatory symmetry between the upper and 
lower states imply the presence of a correspond¬ 
ing phase synchronization between pairs of (a, a 1 ) 
states with the Hanle type factor of Eq. m- 

For our case (1) we therefore set the propor¬ 
tionality constant in Eq. G2D to its maximum 
value that represents full phase synchronization 
whenever m a ' = m a , zero otherwise. This choice 
brings excellent agreement between theory and ob¬ 
servations, as illustrated in Figs. |T] and [3] If we 
would set the proportionality constant to its max¬ 
imum value for all (a, a') combinations we would 
get large positive Q polarization, in contradiction 
with the laboratory experiment, while omitting 
the ground-state coherences leads to exactly zero 
for all wavelengths in Fig. |I] 

Now let us consider our case (2), when the in- 
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cident radiation is 100 % circularly polarized and 
the magnetic field is transverse. In this case the 
electric field of the incident wave packets oscillates 
with equal amplitudes parallel and perpendicular 
to the magnetic field, which means that both the 
q = 0 and q = ±1 resonators of the atomic system 
are excited. This leads to radiative coupling be¬ 
tween all the resonators and therefore between all 
the (a, a') combinations. Consequently we should 
set the constant of proportionality in Eq. (fl^l) to 
its maximum value for all these combinations. Do¬ 
ing this we indeed get excellent agreement between 
theory and experimental data, as illustrated in 
Fig. SI while other choices do not lead to com¬ 
parable success. If we remove all ground-state co¬ 
herences, we get an order of magnitude too small 
polarization. This residual polarization is exclu¬ 
sively due to quantum interferences in the excited 
state. 


3.4. Phase closure as selection rule 


The ground-state coherences that govern the 
polarization of the scattered radiation via Eq. © 
have contributions from both the initial and final 
states. In the previous section we have indicated 
how the initial-state coherences are determined by 
the radiative coupling between the initial m states. 
In the present subsection we will show that all 
final-state interferences are not allowed but that 
they are governed by selection rules that couple 
them to certain combinations of initial-state co¬ 
herences. Therefore an initial-state coherence will 
only contribute to the polarization of the scattered 
radiation if it is coupled in the scattering process 
to a final-state coherence in a way that obeys the 
selection rule. 


The selection rule that governs which combina¬ 
tions of initial and final m states are allowed can 
be formulated as a requirement of ph ase closure , 
a concept that was introduced in Istenflo ( 20091) . 
It has its origin in the geometrical £ factors in 
Eq. ©, which represent the scalar product be¬ 
tween the transverse radiation-field basis vectors 
e Qj p and the complex spherical basis vectors e q 
(with q = 0,±1) of the atomic system. 

The basis vectors e g are defined such that e 0 
is along the magnetic field or quantization axis. 
In this coordinate system the scattering geometry 
is specified by colatitude 9 and azimuth <j>, defined 
with respect to some arbitrarily chosen zero point. 


The scattering geometry of our laboratory exper¬ 
iment is such that 9 = 90° for both the incident 
and scattered beam, while the scattering angle is 
given by <j>' — </>, the azimuth difference between the 
scattered and incident radiation. All the azimuth 
dependence is contained in £^. and e± in the form 
of a phase factor exp(±i</>') or exp(±i<^), re gardless 
of th e values of a and /3 (cf. Eq. (3.86) in IStenflol 


199 4). The complex value of this phase factor de¬ 


pends on the arbitrary choice of zero point for the 
azimuths in the coordinate system used. 


It would be unphysical if the polarization of the 
scattered radiation would in any way depend on 
this choice of coordinate system. Therefore only 
combinations of £ factors for which the choice of 
zero point cancels out can represent physical scat¬ 
tering processes. This gives us a selection rule that 
we refer to as the requirement of phase closure. 
The zero point of the azimuth scale cancels out if 
all the physical effects depend exclusively on the 
difference <f>' — <f> between the scattered and inci¬ 
dent azimuth directions, but never on the absolute 
value of either azimuth. 


Each term in the Mueller scattering matrix is 
made up of the bilinear products between scatter¬ 
ing amplitudes because of the tensor product in 
Eq. ©. Since each scattering amplitude contains 
two £ factors, since e@* = ~^L q when q = ± 1 , 
and since q represents the difference between the 
lower and upper m quantum numbers of the tran¬ 
sition, each Mueller matrix element contains four 
£ factors of the type 


8 OL OL 8 

ff r r ff 

m a -m(, ° mb—mf m^/ —m b , ^m b / —m a r 


(14) 


£ factors relating to the a —> b transition contain 
the azimuth cj> for the incident radiation, while £ 
factors relating to the b —> f transition contain 
the azimuth <f>' for the scattered radiation. If the 
combined complex phase factor for our product of 
four £ factors is exp(*<5), then we see from Eq. (fT4l) 
that 


6 = (m a - mb) 4> + ( m b - mf) </>' 

+(m// - my) 4>' + {jn b ' - m a >) <f >. (15) 

If we would change the zero point and replace 
<j> and <jf with (f> — 4 >o and <jf — tpo , we see that the 
dependence of the phase on <f >o only disappears if 
the factor in front of </> equals the negative of the 
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factor in front of <j>'. This gives us the selection 
rule 

TO/ — TO/' = Too — TO-o' ■ (16) 

This phase closure relation thus demands that 
only the final-state coherences that have the same 
A to as the initial-state coherences are allowed. 
Note that the magnetic quantum numbers of the 
excited state are not involved at all, the rule only 
dictates how the initial and final state coherences 
are coupled. 

In our efforts to reproduce the observational 
data we discovered an additional selection rule 
that seems to govern the coupling between the ini¬ 
tial and final state coherences: 

F f - F f , =F a - F a r . (17) 

It was not anticipated and only found by trial and 
error: With it we got excellent fits to the observa¬ 
tional data, without it the disagreement with the 
data was very large. With afterthought Eq. (ED is 
sensible, since then the scattering transition radia- 
tively couples only coherences that oscillate in the 
same sense. A coherence between F = 1 and 2, 
for instance, oscillates like exp[— i(E\ — E 2 )t/h\. 
A coherence with the reversed order of the states 
would oscillate in anti-phase. Symmetry of the 
scattering process seems to demand that couplings 
with anti-phase oscillations are prohibited. At 
present we lack a full understanding of this rule, 
but it is empirically required to get agreement with 
the observations. 

4. Theoretical modeling of the experimen¬ 
tal data 

Although the calculation of scattering ampli¬ 
tudes from Eq. ([5]) is well-known and uncontrover- 
sial, it is technically difficult for our Di - D 2 case, 
because we have hyperfine structure splitting and 
are in the Paschen-Back regime, in which the Zee- 
man splittings and the transition amplitudes vary 
with field strength in a non-linear way. In support 
of the present project S.V. Berdyugina and D.M. 
Fluri (priv. comm.) kindly provided us with an 
IDL routine that they had written, which calcu¬ 
lates the set of transition amplitudes, resonance 
wavelengths, and energy levels for any given field 
strength in the case of the potassium or sodium Di 
or D 2 system. The correctness of the program has 
been tested through comparison with other codes 


as part of the work of Sowmva et al. (|2014f) . With 
this program the calculation of scattering ampli¬ 
tudes according to Eq. becomes a straightfor¬ 
ward task. 


When trying to theoretically reproduce the ob¬ 
servational data from our laboratory experiment 
there are three other effects to be considered, 
which we have not discussed so far: (1) Effects 
of collisions on the scattering polarization (will be 
dealt with in Sect. 14.11) . (2) Optical depth effects 
with magnetically induced dichroism (will be dealt 
with in Sect. 14.21) . (3) Optical pumping of the 

lower to state populations. 

While (1) and (2) will be accounted for, we 
choose to leave (3) aside in our modeling presenta¬ 
tion here. The reason is that optical pumping that 
creates unbalanced to state populations only mod¬ 
ifies the results quantitatively by modest amounts, 
not qualitatively, in contrast to the fundamental 
issues that we have addressed in the previous sec¬ 
tions concerning the treatment of the initial and 
final state coherences and the selection rules, all 
of which have dramatic effects on the modeling. 
Since the ‘waiting time’ between successive radia¬ 
tive absorption events is in general much longer 
than the time between collisions with the argon 
buffer gas, much of the population imbalances can 
be expected to be erased before each scattering 
event. 


We nevertheless did modeling experiments also 
with inclusion of optical pumping that gives us the 
populations | c a | 2 , but even when the effects of col¬ 
lisions were entirely disregarded in the pumping 
problem the effects on the scattering polarization 
were not dramatic. It was easier to find a good 
fit to the observations when such pumping was 
not included. Pumping does not play an essen¬ 
tial role in the context of identifying the missing 
physics in polarized scattering. Including it intro¬ 
duces new free parameters because the effect of 
collisions on the pumping process, which certainly 
plays a significant role, is not sufficiently known. 
For such reasons attempts to properly deal with 
optical pumping would distract from the main is¬ 
sues of the paper and not lead to improved fits of 
the data. 
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4.1. Effects of collisions 


The vapor cell used for our scattering experi¬ 
ment is small and very convenient to handle. A 
heating coil in the stem of the cell vaporizes the 
metallic potassium to produce the gas in the cell 
center at which the scattering takes place. With 
a gas temperature of 100° C we get high enough 
number density of potassium atoms to provide the 
scattered light needed to measure small degrees of 
polarization. 

To suppress diffusion of the potassium gas, 
which could lead to deposits on the cooler en¬ 
trance and exit windows and make them opaque, 
an argon buffer gas is used in the cell. Unfortu¬ 
nately this has the disadvantage that the potas¬ 
sium atoms are subject to a high collision rate, 
which affects the scattering polarization in signif¬ 
icant ways. 

The collisions cause line broadening by chang¬ 
ing the effective damping constant of the profile 
functions, as well as depolarization of the emit¬ 
ted radiation. Let us, as is customary, denote the 
collisional damping constant by 7 C , while the nat¬ 
ural, radiative damping constant is 7 jv- Then the 
effective 7 that governs the line shape in Eq. m 
is 7 at + 7c- 

The collisional depolarization effects, which de¬ 
stroy the phase relations in the scattering pro¬ 
cess to make the scattered light unpolarized and 
isotropic can be expressed in terms of the depolar¬ 
ization factor 


k c 


In 

7jv + DW ’ 


(18) 


where with AT = 1 or 2 represents the de¬ 

struction rates of the 2K multipoles of the atomic 
system. In ti me-depend ent cla ssical scattering 
theory ( Bommier fc Stenflol 1999 1 


D (K] = 0.5 7 c 


(19) 


for both K = 1 (atomic orientation) and K = 2 
(atomic alignment). In the classical theory the col¬ 
lisions truncate the exponentially damped atomic 
oscillations. After Fourier transformation we then 
get the factor of 0.5 in Eq. (fTf)l) . In detailed quan¬ 
tum mechanical calculations of the scattering pro¬ 
cess the resulting value for this factor depends on 
the atomic system but stays typically within at 
least a factor of two around the classical value. In 


the absence of other data our best choice is there¬ 
fore to use the relation represented by Eq. (USD. 

The depolarization factor k c is used to scale 
all elements of the Mueller scattering matrix ex¬ 
cept the first one (Mu), which represents isotropic 
scattering. This is a simplified treatment, but in 
view of the large collisional rates and the use of 
these rates for a parametrization of the problem, 
it is not justified here to go into the intricacies of 
partial frequency redistribution theory. 

The suppression effect from the depolarization 
caused by k c is largely compensated for by the 
increased overlap and interference between the 
ground-state sublevels due to the collisional broad¬ 
ening, as described by Eq. m The question 
is however what collision-modified expression for 
7 we should use in Eq. m, e.g. 7 N + 7 C or 
7 at + ? On the one hand one may argue for 

the first of these two choices, since these coher¬ 
ences are determined by the radiative couplings 
between the various combinations of allowed tran¬ 
sitions, and these couplings are governed by the 
profile functions in which 7 = yjv + 7 C - On the 
other hand there is a good case for the second 
expression, since the depolarization factor used 
for the upper-state coherences is based on D^ K \ 
Through experimentation with both versions we 
find that the best fit results are obtained if we 
use an average between the two cases. This im¬ 
plies that we for the effective 7 in Eq. m use 
7 jv + 0 . 757 c . For symmetry reasons we use the 
same expression for the initial and final state co¬ 
herences. 

Let us stress that the choice of the expression 
for 7 to be used for the ground-state coherences 
is not crucial for our investigation. It does not 
lead to any qualitative changes, but only affects 
the quantitative details of the model fitting. Our 
aim here is to make the choices that appear most 
reasonable while minimizing the number of free 
parameters. 

If we thus fix the choices for the three re¬ 
lations used to describe the effects of collisions, 
7jv + 7c to describe the shape of the profile func¬ 
tion, 7 jv + 0 . 57 c to describe the collisional depo¬ 
larization factor used for the Mueller matrix ele¬ 
ments, and 7 n + 0 . 757 c to describe the effect on 
the ground-state coherences, then we are left with 
a single free collisional parameter, namely 7 C . 
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4.2. Optical depth effects 

The optical depth of the potassium gas must 
be non-zero, otherwise we would not get any scat¬ 
tered photons. Choosing the working temperature 
of the cell represents a compromise: we want to 
optimize the number of scattered photons to get a 
good signal, while avoiding multiple scattering or 
excessive optical depth effects. The finite optical 
depth leads to enhanced absorption at line cen¬ 
ter, which may change the profile shape, and also 
causes dichroism when a magnetic field is present. 
Comparison between our theory and the experi¬ 
mental data shows that the optical depth effects 
are too small for significant multiple scattering, 
but large enough to cause significant line broad¬ 
ening and magnetically induced dichroism. They 
therefore need to be accounted for, but a single 
free parameter, the optical depth r at line cen¬ 
ter (where the optical depth reaches its maximum 
value) is sufficient. We find that the optimum fit 
value for r is 0.25 in the case of Di, while it is 
twice this value for D 2 because of its twice higher 
oscillator strength. 

Our scattering experiment can be described as 
a transformation of the input Stokes vector Si„ 
into the output Stokes vector S out through 

S ou t = M e ff S in , (20) 

where M c ff is the effective Mueller matrix of the 
system. When optical depth effects are accounted 
for, it can be represented by 

Meff =M arm M scat M 

arm ? (21) 

where M scat is the scattering Mueller matrix, 

the same as the M that is given by Eq. ©. 

while M arm represents the Mueller matrix for the 
medium in an arm of the cell. There are two rel¬ 
evant arms here: The incident light goes through 
the input cell arm before it reaches the cell center, 
where it is scattered and then traverses the output 
cell arm before being analyzed for its polarization 
state and being detected. The Mueller matrices 
for the input and output arms are the same and 
surround the scattering matrix. 

In the unpolarized case with an optical depth 
r at line center and an absorption profile <f> that 
is normalized to unity at line center, the effect 
of the medium would be absorption by the factor 
exp(—</>r) (since the contributions from emission 


through scattering in the pencil-shaped geometry 
of the medium would give a negligeable contribu¬ 
tion in the particular beam direction in compari¬ 
son with the absorption effects in that direction). 
In the polarized case, when the medium in a cell 
arm is described by the 4x4 Mueller matrix $ 
instead of by a scalar absorption coefficient, the 
dichroic effects of the medium are represented by 

M arm =e-* r . (22) 


The exponentiation of a matrix is defined by its 
Taylor expansion. 

In a magnetized medium 


< <t>I 

0Q 

4>u 

<t>v\ 

<t>Q 

4>i 

Vv 

-ipu 
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-Vv 
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where 4>q,u describe the transverse Zeeman effect, 
(j>y the longitudinal Zeeman effect, and 4'Q,u,v 
the magneto-optical effects, while </>/ represents 
the effect of the Zeeman splitting on the inten- 
sit y absor ption profile (cf. Eqs. (6.58) and (6.59) 


Stenflo 


ption p 


Since we are dealing with a very weakly magne¬ 
tized medium (all field strengths that we consider 
here are below 20 G), all off-diagonal elements of 
$ are small in comparison with </>/. In our power 
series expansion of exp(—$ r) we therefore only 
need to retain the terms of first order in these ele¬ 
ments, while allowing the optical depth r to be of 
any higher order. 


Here we will only model the case when the in¬ 
cident radiation is propagating perpendicular to 
the magnetic field direction, which also defines the 
positive Stokes Q direction. Then all dichroic ele¬ 
ments of are zero except </>/ and (j)Q. In the case 
of circularly polarized light plays no role, the 
whole effect is in the form exp(—</>/ r). However, in 
the case of linear Q polarization, (f>Q is significant 
and represents the effects of magnetically induced 
dichroism. 


Without this dichroism we are unable to ex¬ 
plain the observed dependence of the scattered Q 
polarization on field strength, since our scattering 
theory alone does not predict any significant ef¬ 
fect, even when the ground-state coherences are 
included. Comparison with our data gives us an 
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estimate of the value of r that is needed to explain 
the observed field dependence. Another constraint 
comes from the observed profile shapes. The line 
width increases with increasing r. It turns out 
that the value of r that gives us the right field de¬ 
pendence automatically also gives us the right line 
width. 

For the calculation of the 4> matrix we have 
to sum over the contributions from the various (f> q 
profile functions (q = 0, ±1) weighted by the corre¬ 
sponding transition strengths t 2 abl to arrive at the 
expressions for 4>j and <pQ. In this case the profile 
functions (f> q are described by Voigt functions, hav¬ 
ing a Doppler width given by the thermal width at 
100 ° C (10.2 mA) and a damping constant given 
by 7 jv + 7 c i which is the same damping rate that 
governs the profile functions of the scattering pro¬ 
cess. 

4.3. Using D 2 to determine the free pa¬ 
rameters 

With our previous definitions we are left with 
only two free model parameters, 7 C , which governs 
all the collisional effects, and r, which governs the 
diclrroic or optical-depth effects. For physical con¬ 
sistency we further have the constraint that the 
same parameter values must be used for both the 
Di and D 2 observations. We are not allowed to 
use independent fit parameters for these two lines, 
since the same vapor cell is used for both. The 
only modification is that we need to use 2r for 
D 2 , since it has twice as large oscillator strength 
and therefore twice as large optical thickness. 

The ground-state coherence effects that we have 
been discussing play no significant role for the D 2 
line, because D 2 has intrinsic polarization (from 
the upper level) that is larger by approximately a 
factor of 50 and therefore masks the ground-state 
effects. For this reason the D 2 line can perfectly 
well be described by standard scattering theory. 
By using the D 2 line to fix the values of our two 
free parameters, the obtained values become inde¬ 
pendent of the way in which we are dealing with 
the ground-state coherences. No further parame¬ 
ter fitting is then needed when addressing the Di 
problem. 

Figure [2] shows the excellent fit we get to the 
observed Stokes I and Q D 2 profiles for the non¬ 
magnetic case when the incident radiation is 100 % 


linearly polarized perpendicular to the scattering 
plane (Q type polarization). The dashed curves 
represent the observations, the solid curves the 
theoretical model with y c = 90 jn ~ 3.5 x 10 9 s -1 
and r = 0.25 (implying that the maximum optical 
depth at the D 2 line center is 0.50). The theo¬ 
retical profiles have been convolved with a Gaus¬ 
sian having a Doppler width that equals the ther¬ 
mal width for 100° C. No additional instrumental 
broadening has been applied (we have no reason 
to expect any significant spectral smearing of the 
laser beam). 

The fit in Fig. [2] is nearly perfect, except that 
there are some deviations in the distant wings of 
the I profile. The Q polarization amplitude scales 
directly with the depolarization factor k c , which 
therefore gets fixed precisely by the data. Through 
Eqs. (fTHl) and (fT9l) this determines the value of 7 C . 

The line width scales with j c , but it also in¬ 
creases with increasing r, since the absorption 
mainly suppresses the line core while having lit¬ 
tle effect on the wings. If we use the y c that is 
fixed by the Q amplitude but set r = 0, then the 
theoretical profiles are significantly narrower than 
the observed ones, and there would be a need for 
additional line broadening (e.g. due to some hypo¬ 
thetical instrumental smearing) to get agreement 
between the profile shapes. With our chosen value 
for r we get an excellent profile fit without the 
need for any such additional ad hoc broadening. 

As we will see in the next subsection, the r that 
is needed to fit the D 2 profile shape is also the r 
that is needed to account for the field strength de¬ 
pendence of the Di Q polarization. It is gratifying 
to note that there is such consistency between the 
independent D 2 and Di constraints on r. 

If we would increase the value of r significantly 
further, so that 2r would start to become compa¬ 
rable to unity, then the D 2 Stokes I profile not only 
broadens but gets deformed by becoming double¬ 
humped with a depression at line center. The ab¬ 
sence of any sign of such a deformation shows that 
we can safely rule out the possibility that the op¬ 
tical depth could be much larger than our present 
estimates. This conclusion refers only to the quali¬ 
tative profile shape without using numerical values 
of any amplitudes or widths. 
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4.4. Role of magnetically induced dichro- 
ism 

In Fig. [Tjwe showed the results for Di scatter¬ 
ing in the non-magnetic case with the same ge¬ 
ometry and same parameter values as for the D 2 
plot in Fig. [2] We note that the D x polarization is 
of opposite sign and with about 50 times smaller 
amplitude than that of D 2 . The laboratory exper¬ 
iment has the sensitivity to measure polarization 
effects down in the 10 -5 range. Such a sensitivity 
is needed to uncover the tiny polarization effects 
and profile shapes for Di scattering. The access 
to this parameter range allows us to explore scat¬ 
tering physics in previously untested domains. 

When for Di doing the Stokes Q scattering for a 
sequence of different field strengths (with the field 
oriented perpendicular to the scattering plane) the 
laboratory experiment reveals a significant field 
dependence that scales with the square of the field 
strength. In contrast, a significant field depen¬ 
dence (within the studied field strength range) 
is absent in the computed Q polarization when 
the effects of the surrounding medium are not ac¬ 
counted for. If we however introduce a finite opti¬ 
cal depth, then we get a field dependence that also 
scales with the square of the field strength, as it 
should when it is due to the transverse Zeeman ef¬ 
fect. The amplitude of this dependence scales with 
r. Choosing r = 0.25 gives good agreement with 
the observational data, as illustrated by Fig. [3] 

The main source of this field dependence comes 
from the <f>Q element of the dichroic 4> matrix of 
Eq. (l23ll . which represents the transverse Zeeman 
effect. In absorption the contribution from the it 
component in the line core is negative, while the 
a contributions in the wings are positive, but in 
the <f>Q profile they are balanced to make the inte¬ 
grated profile zero. However, the dichroic contri¬ 
bution to the Q polarization comes not from cj>Q 
alone, but from its product with the intensity pro¬ 
file. This gives much more weight to the negative 
contribution from the line core than to the wing 
contributions. Therefore the net contribution is 
negative, with the consequence that the already 
negative non-magnetic scattering Q polarization 
gets enhanced. 

Figure [3] shows the I and Q profiles for the 
largest field strength used in our laboratory ex¬ 
periment, 18.7 G. As in Fig. |T] the observations 
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Fig. 2.— Stokes I and Q profiles for potassium 
D 2 in the non-magnetic case when the incident 
radiation is 100% linearly polarized in the Q di¬ 
rection (perpendicular to the scattering plane) like 
in Fig. [I] that represented the Di case. The same 
values for the free parameters, q c = flOqjv and 
r = 0.25, are used for both the D 2 and Di cases. 
Note that the D 2 polarization is about 50 times 
larger than the Di polarization and of opposite 
sign. 
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Fig. 3.— Stokes I and Q profiles for K Di for the 
same geometry and parameters as in Fig.[T] except 
that here we show the case when the vertical mag¬ 
netic field is 18.7 G (while it was zero in Fig. [TJ. 
Dashed curves: observations; solid curves: theory. 
The difference in polarization scale with respect 
to Fig. |T] is due to magnetically induced dichroism 
in the potassium gas. The systematic wavelength 
shift of the dashed Q curve relative to the solid 
one can be understood in terms of the time con¬ 
stant used for the lock-in amplifier to optimize the 
polarimetric sensitivity. 
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are represented by the dashed curves, the the¬ 
ory by the solid curves. The geometry and other 
parameters are the same as for Fig. [lj the only 
difference being the field strength. We see that 
an excellent fit is obtained, which would not have 
been the case without introducing the finite opti¬ 
cal depth r = 0.25. Comparison between Figs. [3] 
and Q] shows that the profile shapes are almost un¬ 
affected by the magnetic field, but the polarization 
scales in the two figures are significantly different 
(although the difference is rather modest). 

Both Figs. [3] and Q] give the (false) impression 
that the theoretical model is unable to reproduce 
the exact wavelength position of the Q profiles. 
The observational (dashed) curves for the Q pro¬ 
files are seen to be systematically shifted towards 
the right relative to the theoretical (solid) curves. 
The magnitude of the shift is about 2.8 mA, which 
corresponds to about 1.7 step sizes in the stepwise 
wavelength tuning of the laser to scan the line pro¬ 
file. 

The probable origin of this shift has to do with 
the time constant used by the lock-in amplifier 
when demodulating the polarization information, 
which for Stokes Q is at 84 kHz due to the piezo- 
clastic modulation. To achieve polarimetric sensi¬ 
tivity down into the 10~ 5 range one needs maxi¬ 
mum temporal integration of the demodulated AC 
signal before it is read out. The chosen analog in¬ 
tegration time should however be compatible with 
the scan rate. If the integration time is longer 
than the time between tuning steps, there will be 
a delay and overlap between the steps, resulting in 
shift and broadening of the demodulated profile. 
As such an analog time constant is not applied 
to the DC (Stokes I) signal, the Q profile will be 
shifted (and slightly broadened) in the tuning di¬ 
rection with respect to the / profile. The observed 
shift that is of order 1-2 tuning steps is a natural 
consequence of trying to maximize the polarimet¬ 
ric sensitivity of the measurements. 

Note that the time constant could be much re¬ 
duced in the D 2 case, since the polarization signals 
were about 50 times larger. Therefore we do not 
see any corresponding shift of the Q profile in the 
plot for D 2 in Fig. [2] 


4.5. Circularly polarized scattering for 
transverse magnetic fields 

Let us now turn to the Di case when 100% 
right-handed circularly polarized incident radia¬ 
tion is used, and Stokes I and V are measured 
in the scattered radiation. The geometry remains 
the same as in our previous discussion of Q scat¬ 
tering, with the magnetic field being perpendicular 
to the scattering plane. 

The laboratory measurements were made for a 
sequence of field strengths between ±18.7G. In 
Fig. []] we represent the results as the fractional 
circular polarization V/I (%) vs. the strength of 
the vertically oriented field. In this case V rep¬ 
resents the wavelength integrated V, I the inte¬ 
grated Stokes I profile. The observations are given 
by the filled circles and the dashed analytical fit 
function, while the solid line is obtained with our 
theory, using the same values for parameters q c 
and t as we used for the modeling of Q scattering. 
In contrast to the linearly polarized case, the in¬ 
cident circular polarization leads to radiative cou¬ 
plings between all the initial m states and not only 
between those with identical m quantum numbers, 
as explained in Sect. 13.3.31 

We notice in Fig. [4] the excellent agreement be¬ 
tween the solid and dashed curves, in particular 
the equal slopes of the field dependence, except 
for fields in the inner range 0-2 G. The remark¬ 
able negative peak is centered at the value 0.43 G 
(marked by the vertical dotted line in the dia¬ 
gram). This offset from zero can be understood as 
being caused by the Earth’s magnetic field, which 
has a downwards pointing component of approxi¬ 
mately this magnitude. To compensate this com¬ 
ponent our imposed magnetic field needs to have 
an upwards (positive component) of this magni¬ 
tude, which leads to the apparent small shift of 
our field strength scale. 

Since the Earth’s magnetic field is not exactly 
vertical but inclined, the residual field will become 
horizontal when the vertical field gets compen¬ 
sated, but the azimuth orientation of our exper¬ 
imental setup with respect to the azimuth of this 
residual horizontal component is not sufficiently 
known to allow meaningful modeling of this be¬ 
havior. In the milligauss range in the immedi¬ 
ate neighborhood of the +0.43 G offset we can ex¬ 
pect the magnetic-field vector to rapidly change its 
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Fig. 4.— Degree of circular polarization in the 
beam scattered at 90° by potassium Di, plotted as 
a function of field strength. The fractional circular 
polarization here represents the area of the Stokes 
V profile divided by the area of the / profile. The 
incident beam has 100 % right-handed circular po¬ 
larization, and the magnetic field is perpendicular 
to the scattering plane. The observations are rep¬ 
resented by the filled circles and the dashed curve, 
the theoretical model that includes initial and fi¬ 
nal state coherences by the solid line. When these 
ground-state coherences are omitted but the inter¬ 
ference terms of the excited state are retained, one 
obtains the dotted line, which has a slope that is 
an order of magnitude smaller. If all m state co¬ 
herences would be omitted the polarization would 
be exactly zero. 


orientation between vertical and horizontal. The 
sharp negative peak in Fig. 0] shows that the scat¬ 
tering polarization has enormous sensitivity to the 
field orientation in this milligauss range. The ex¬ 
istence of this sharp peak indicates that our de¬ 
scription of the physical system is still incomplete. 
The radiation-matter interaction appears to have 
a much richer substructure than anticipated. 


Other laboratory tests that we have made show 
that when we have horizontal magnetic fields 
aligned along the propagation direction of the 
incident radiation and choose this radiation to be 
100 % circularly polarized, then not only the scat¬ 
tered Stokes V, but also the scattered Stokes I, 
varies dramatically with the strength of the im¬ 
posed field. In contrast to the vertical field case, 
for which Stokes / is insensitive to field strength, 
the horizontal field case leads to a Stokes / that 
has a sharp peak around zero Held strength and 
rapidly de clines towards zero as t he field strength 
increases ( Thalmann et al. 20091) . This behavior 


can be readily understood in terms of extreme op¬ 
tical pumping of the medium. With the described 
geometry either only q = +1 or q = —1 transi¬ 
tions will be excited (depending on the handed¬ 
ness of the incident circular polarization). If for 
example only radiative excitations that raise the 
m quantum number by +1 are possible, then all 
the atoms of the medium will quickly be pumped 
into the ground state with m = +2, F = 2, from 
which no more radiative excitation is possible. 
This means that the medium becomes transpar¬ 
ent and the scattered intensity goes to zero, as 
observed except for the weak fields around zero. 
The Stokes I peak at zero field strength occurs 
because there are weak stray fields (in particular 
the Earth’s magnetic field), which break the sym¬ 
metry and prevent the extreme pumping scenario 
from happening. 


The laboratory experiment was also done with 
the opposite handedness of the input polarization 
(100% left-handed circularly polarized radiation). 
The result is the same, except that the sign of 
the output V/I is reversed (and consequently the 
sharp peak at +0.43 G is positive). The exact sign 
reversal also happens for the theoretical model, 
the diagram looks the same if we reflect it in the 
horizontal axis. This symmetry behavior implies 
(as verified by our theory) that all the V signal 
in the output radiation comes from the non-zero 
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Mueller matrix element M 44 , while the other ele¬ 
ments of the 4th row, like M 41 , are exactly zero. 

In contrast to the scattering of Q polarization, 
dicliroism does not play any role in Fig. [I] The 
reason is that the term <\>y , which represents the 
longitudinal Zeeman effect in the Mueller matrix 
is zero because we have chosen the magnetic 
field to be transverse. Therefore the effect of the fi¬ 
nite optical depth is to multiply both the scattered 
I and V profiles with the same factor, exp(—</>/ r), 
which differentially suppresses the line core more 
than the wings. Since it turns out that the scat¬ 
tered / and V profiles have identical shapes and 
thus are deformed in the same way by the factor 
exp (—(f>i t), the fractional polarization, integrated 
V divided by integrated /, is not affected. This is 
not so in the Q scattering case, since the I and 
Q profiles have different shapes, and in addition 
there is the contribution from 0 q, which is non¬ 
zero since it represents the transverse Zeeman ef¬ 
fect. 

If we would remove the initial and final state co¬ 
herences but keep the upper-state coherences (like 
in the standard scattering theory of Eq. ©), then 
we obtain the dotted curve in Fig. [H which has a 
slope that is an order of magnitude smaller than 
that of the observations. If we would also remove 
the upper-state interferences we would get exactly 
zero for all field strengths. This demonstrates 
that all the polarization effects in V scattering 
for transverse magnetic fields have their origin in 
coherences between separate magnetic substates, 
and that the dominating contributions come from 
coherences in the ground state. 

5. Limitations, conclusions, and outlook 

Our search for missing physics in quantum scat¬ 
tering was motivated by the observed existence of 
a small but significant polarization peak in the 
Sun’s spectrum of the Nai Di line at 5896 A. 
Quantum theory appeared to predict that the Di 
line should be intrinsically unpolarizable, and the 
transverse Zeeman effect could be ruled out as a 
source of the observed polarization. After a decade 
of failed attempts to explain this enigmatic feature 
we decided to set up a laboratory experiment to 
address the question whether this was a problem 
of solar physics or quantum physics. 

Quantum systems often contain an inherently 


rich and complex structure that is not always 
fully understood. The correct application of a 
complex theory needs guidance from experimen¬ 
tal data. The problem is that the theory of polar¬ 
ized scattering at multi-level atomic systems has 
never been experimentally tested before in the pa¬ 
rameter range that we are addressing here. This 
needs to be done before we can have sufficient con¬ 
fidence in the application of the theory to ‘messy’ 
astropliysical environments over which we have no 
control. 

In most cases when new types of instrumenta¬ 
tion allows the exploration of parameter domains 
that have never been accessed before, surprising 
new phenomena have been discovered which were 
not predicted by the theory available at the time. 
The Di enigma was never expected but was a 
major surprise when the Sun’s spectrum could 
be recorded with high-precision imaging spectro- 
polarimetry allowing a sensitivity of 10~ 5 in the 
fractional polarization. Laboratory experiments 
in polarized scattering was a hot topic in the 
early days of quantum mechanics but they were 
largely abandoned around 1935 in favor of other 
topics. The polarization phenomena that we are 
discussing here are way out of reach for the exper¬ 
imental technology available at that time. 

The aim of our laboratory experiment is not to 
try to emulate scattering in the solar atmosphere, 
but to test the basics of Di scattering physics un¬ 
der simple and controlled conditions. The exper¬ 
iment revealed that there is indeed a problem of 
quantum physics: Di exhibits a rich polarization 
structure in situations where available theory pre¬ 
dicts zero polarization. 

When abandoning the prevailing ‘flat-spectrum 
approximation’ we open a Pandora’s box to a new 
world of coherency phenomena that were excluded 
before. In the present paper we show that polar¬ 
ized scattering theory for multi-level atomic sys¬ 
tems has to be extended to include both the initial- 
state and final-state coherences (in addition to the 
coherences in the excited, intermediate state) in 
order to explain the observed polarization struc¬ 
tures. With our phenomenological extension of 
the theory we succeed in reproducing the obser¬ 
vations not only qualitatively in terms of profile 
shapes and signs, but also quantitatively in great 
detail, both for the scattering of linear and circular 
polarization in a transverse field geometry. This 
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gives us confidence that the missing ingredients of 
the theory have been identified correctly. 

There are however also indications that our the¬ 
oretical extension is far from complete. Experi¬ 
ments with the scattering of circular polarization 
indicate that there are dramatic effects happen¬ 
ing in the milligauss regime of the external mag¬ 
netic field (cf. Fig. Q}. We are presently unable to 
deal with this regime, because we lack control and 
knowledge of the strength and orientation of the 
field in the milligauss range due to the presence of 
stray fields, including the Earth’s magnetic field. 
An experiment that could properly deal with this 
issue is likely to reveal additional substructures of 
the atomic system that have not been identified so 
far. To fully address the behavior for such weak 
fields the experiment needs to be set up in a mag¬ 
netically clean environment, which is not easy but 
doable. 


Another deficiency is that we only had access 
to a vapor cell that used argon buffer gas to pre¬ 
vent condensations on the cooler cell windows. As 
a consequence the potassium atoms experience a 
high collision rate (with 7 C « 90 7 jv), which pro¬ 
duce large line broadening and depolarization ef¬ 
fects. These effects had to be accounted for in a 
phenomenological and parametrized way. In fu¬ 
ture one would like to do the experiment with a 
cell that does not use any buffer gas, so that one 
can explore the clean case of scattering in the colli¬ 
sionless domain. In the solar chromosphere, where 
the core of the Nai Di line is formed, one should 
approach nearly collisionle ss co nditions, but ac¬ 
cording to Kerkeni fc Bommier (2002fl lower level 
polarization in Nai Di will always be destroyed 
under solar conditions. Ideally one would like an 
experimental setup that would allow exploration 
of the scattering polarization as a function of col¬ 
lision rate, by being able to change and control the 
amount of buffer gas being used. 


It is generally problematic to apply an estab¬ 
lished physical theory to parameter domains in 
which the theory has never been tested by experi¬ 
ment. We do not have a theory that has been suf¬ 
ficiently tested experimentally when it comes to 
the interpretation of scattering polarization from 
multi-level atomic systems. Experimental explo¬ 
ration of this poorly studied domain may bring 
unexpected new insights into neglected aspects of 
quantum physics. In the present paper we have 


tried to indicate the direction in which the explo¬ 
ration of this domain could take us. 

The potassium cell was specially tailor made 
free of charge for our laboratory experiment by 
the late Alessandro Cacciani, who had devoted 
much of his life to perfect the art of making such 
cells for magneto-optical filter systems used in he¬ 
lioseismology and solar magnetometry. The scat¬ 
tering experiment was set up and carried out at 
ETH Zurich as part of two PhD projects, by Alex 
Feller and Christian Thalmann. The data bank 
set up by Thalmann for the collection of his labo¬ 
ratory recordings has served as a resource to guide 
the theoretical efforts. I am grateful to Svetlana 
Berdyugina and Dominique Fluri for providing 
me with their IDL code to calculate the atomic 
level structure and transition amplitudes in the 
Paschen-Back regime for the sodium and potas¬ 
sium Di - D 2 systems, and to Roberto Casini for 
various discussions on quantum physics, in partic¬ 
ular on the nature of the flat-spectrum approxi¬ 
mation. 
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